Achieving a BCS transition in an atomic Fermi gas 
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We consider a gas of cold fermionic atoms having two spin components with interactions char- 
acterized by their s-wave scattering length a. At positive scattering length the atoms form weakly 
bound bosonic molecules which can be evaporatively cooled to undergo Bose- Einstein condensation, 
whereas at negative scattering length BCS pairing can take place. It is shown that, by adiabatically 
tuning the scattering length a from positive to negative values, one may transform the molecular 
Bose-Einstein condensate into a highly degenerate atomic Fermi gas, with the ratio of temperature 
to Fermi temperature T/Tf ~ 10~ 2 . The corresponding critical final value of fcp|a| which leads to 
the BCS transition is found to be about one half, where &f is the Fermi momentum. 

PACS numbers: 



Much progress has been made in the achievement of 
increasingly degenerate regimes of trapped atomic Fermi 
gases 0,|2SElME10- One °f the major goals of studies 
of these systems is to observe a transition to a paired- 
fermion superfluid state. In a recent MIT experiment on 
a sympathetically cooled single component Fermi gas, a 
system with T/Tp = 0.05 and 3 x 10 7 fermions was real- 
ized 0, where Tp is the Fermi temperature. As trapped 
atomic gases present an essentially impurity free system, 
with density, temperature, and interaction strength all 
free parameters, they offer great opportunities for inves- 
tigation of fundamental theories of superfluid states. 

An important new step in experiments has been the 
use of Feshbach resonances, whereby one may tune the 
s-wave scattering length a from positive to negative over 
many orders of magnitude This opens up the pos- 
sibility of investigating the Bose-Einstein condensate to 
Bardeen-Cooper-Schrieffer (BEC-BCS) crossover |i| 
[H[l3 The BCS limit occurs for fc F |oj < 1, where k F is 
the Fermi momentum and a < 0. In this case the fermion 
pair size is much larger than the interatomic spacing. 
The BEC limit occurs for a positive scattering length 
much smaller than the interparticlc separation. Then 
the fermions form weakly bound dimers of a size ~ a, 
and the Bose-Einstein condensation of these composite 
bosons may be described by the well-developed theory of 
trapped BEC's (l4|. Recently, long-lived molecules have 
been created in a reversible, and therefore adiabatic fash- 
ion from a degenerate Fermi gas by tuning the scattering 
length |l5|. This provides many possibilities for evapo- 
ratively cooling the molecules into the BEC regime. In 
contrast, a BCS transition requires very low tempera- 
tures not yet obtained in two-component Fermi gases. 

In this Letter, we present a straightforward theoretical 
analysis which shows that a deeply degenerate Fermi gas 
may indeed be created from a molecular BEC. By adi- 
abatically changing the scattering length from positive 
to negative values, the entropy is held constant, which 
leads to a strong decrease in temperature. This sug- 



gests a novel method of cooling a Fermi gas to extremely 
low temperatures. One can then reach a temperature 
T ~ f Q~ 2 Tp and achieve a BCS transition. 

Consider a harmonically trapped Fermi gas described 
by the grand canonical ensemble with two equally pop- 
ulated spin states. The entropy of the gas above the 
critical temperature is equal to the entropy of an ideal 
Fermi gas in a trap with mean field corrections. Omitting 
these corrections the grand potential Jl is 01 



dep(e) ln[l - n(e)} 



with p(e) the density of states and 

n(e) = l/{exp[/3(e - //)] 
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the Fermi weighting factor. For a harmonic trap, the 
density of states for a non-interacting two-component gas 
in the semi-classical limit is 



p(e) = e 2 /(hco) 
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Substituting Eq. © into Eq. Q and integrating 
by parts, one obtains the simplified integral SI = 
— (A/3) J °° dee 3 n(e). In the degenerate regime, given by 
jj, 3> fcsT, the integral may be expanded as 
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The entropy may then be obtained from the relation [T^ | 
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giving 



S ~ ks — 



S = -drSl\ 
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A relation between the total number of particles and 
the chemical potential follows from the equation N = 
— dfj,n\T,n and reads 

N~(A/3)[^ + n 2 (k B T) 2 fi}. (7) 
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So, to lowest order, the chemical potential is \i = 
(3N/A) 1 / 3 and the expression for the entropy becomes 



S = k B Nn 2 T/T F + 0(T 3 ) , 



(8) 



where Tp = (3N) 1 / 3 hco is the Fermi temperature of 
the non- interacting gas. In the case of an interacting 
Fermi gas, the lowest order mean field correction in the 
Thomas-Fermi limit leads to the appearance of an ex- 
tra multiple (1 + 64fci?a/357r 2 ) in the numerator of the 
right hand side of Eq. (JHJ. This gives a correction of less 
than 10 percent to the entropy for kp\a\ < 1/2. It may 
therefore be neglected in our calculations. Equation (JSJ 
is valid for k B T > hut in an isotropic potential. Its valid- 
ity extends to much lower temperatures (although much 
larger than the inverse density of states at the Fermi en- 
ergy) for incommensurable trap frequencies. 

Consider the analogous calculation for a BEC of weakly 
bound molecules at a > 0. We assume that the con- 
densate is in the weakly interacting regime, with rj = 
flmoi "C 1 where n mo \ is the molecular density in the 



i 1/3 

l mol 



trap center and the scattering length of two molecules 
a mo i is related to the atomic scattering length a by a mo i = 
0.6a 0|- Note that these weakly bound molecules ex- 
hibit a remarkable collisional stability at large values of 
a fl5Hl8l |. This stability originates from a strong decrease 
in the relaxation rate to deep bound states with increas- 
ing a, and should enable efficient evaporative cooling [lij . 
A small value of r\ may be reached in experiments by re- 
ducing the trap frequencies. 

If r] <C 1, the mean field shift in the condensation 
temperature is <5T BE c/7bec = -3.25T? 5 / 4 Q 

and the 

quantum depletion of the condensate is negligible so that 
Bogoliubov theory may be used. It also implies 



k B T BEC /E B = 7.78T? 5 / 2 < 1 



(9) 



where E B — Ti 2 /ma 2 is the molecular binding energy. 
This condition, together with T < TbeCj ensures that the 
molecules may be treated as weakly interacting bosons 
and any effect due to their thermal dissociation is negli- 
gible. The grand potential is 



n = k B T / dep{e) In (l 
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where f3 now refers to the temperature of the Bose gas. 
The density of states for an interacting BEC in the Bo- 
goliubov approximation may be calculated directly from 
the Bogoliubov Hamiltonian as p(e) = Tr[<5(e — Hbog)]- In 
the Thomas-Fermi limit, which holds for a molecular con- 
densate with chemical potential /x mo i 3> fiuj, the density 
of states may be approximated semi-classically by 



drdp 
(2irh) 3 



S[e - e bog (p,r)} 
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where the energy of Boguliubov excitations is (see Q| 
and references therein) 



_ J yJp 2 /2m mQ \ [p 2 /2 + 2g mo \n Q (r)] \r\ < R 
h ° 9 \ P 2 /2m mo i + |m mo iw 2 r 2 - p mo \ \r\ > R. 

(12) 
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Here m mo i = 2m is the mass of a molecule and 

n a (r) = (/i mo i/.gmoi)(l - r 2 /R 2 ) 

is the Thomas-Fermi density profile, with R = 
(2/Ltmoi/TOmoi^ 2 ) 1 ^ 2 - The coupling constant for the 
molecule-molecule interaction is g mo i — 47r?i 2 a mo i/m mo i. 
The exact calculation of Eq. I|ll|) yields 

P(e) = [MLi/(^ 3 ^ 3 )]{2%/2ztanh- 1 [V2l/(l + z)} 
+4z 3/2 - V2z 2 [n + 2tan~ 1 ((l - z)/V2z)} 



-(1 



9 - sin(40 o )/4]} , (14) 



is the rescaled energy and 9 = 



where z = e/fj, n 
cos _1 (l/Vl + z). In the limit k B T > /i mo i, 
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Given the density of states, the entropy can be calculated 
from Eq. JSJ and (|10|) . The resulting expression is 



S 



) (tbec) 
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where £ is the Riemann zeta function, with £(3) = 
1.202 . . ., N mo \ = N/2 is the number of molecules and 



G(u) 



dzf(z) 



1 



ln(l 



, (17) 



with u = /3/i mo i, and f(z) = [(riw) 3 /^ 2 nol ]p(z^ mo i) so as 
to make the units explicit. Equation i|17|) may be inte- 
grated numerically, or explicitly in the limit of Eq. i|15fl . 
in which case 



/ T \ 3 / 2tt 4 
V Tbec J \45C(3) + 6 ~k^T ) 



(-3^ (18) 

This expansion is accurate to within 10% for k B T / fi mo \ > 
1/10. In the above expression for the entropy, the con- 
densation temperature for a non-interacting Bose gas, 
k B TBEC = ^[Nmoi/^S)] 1 ^ 3 was used. Setting the en- 
tropy of the molecular condensate and of the Fermi gas 
equal one obtains the final temperature of the Fermi gas: 



T_\ = G(/?Mmol) / 



T V 
Tbec / ; 



(19) 
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where (3 refers to the initial temperature. 

In the strongly degenerate regime, evaporative cooling 
of a Fermi gas becomes difficult, so that until recently a 
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FIG. 1: Final temperature of the Fermi gas (left axis) and 
critical condition for a BCS transition (right axis) as a func- 
tion of the initial temperature of the Bose condensed molecu- 
lar gas that has been transformed to a normal attractive Fermi 
gas via adiabatic switching of the scattering length. The 
curves, from upper to lower, correspond to /imoi/fcsTBEC = 1, 
1/2, and 1/4, where jU mo i is the molecular chemical potential. 



maximum degeneracy of T/Tp ~ 0.2 appeared to be a 
lower limit, preventing the observation of the BCS tran- 
sition. The adiabatic tuning of a molecular condensate 
suggests an alternate method to cooling the fermions di- 
rectly. According to Eqs. JHJ) and l|18JI . the entropy is 
proportional to T for fermions and T 3 for bosons. The 
limits of evaporative cooling for a Thomas-Fermi Bose 
condensed gas are given by (i ~ kp,T . Thus the final 
degeneracy that may be obtained by adiabatic switch- 
ing of an evaporativcly cooled molecular BEC is given by 
Eq. with G(j3fM mo {) = G(l) ~ 8.32. Larger values 

of /3/i m oi are in principle achievable, e.g. by sympathetic 
cooling. In Fig. ^ is shown the final temperature as a 
function of the initial temperature for various interaction 
strengths given by the ratio ^ mo i/fcsTBEC- 

The lowest experimentally achieved temperature of a 
degenerate Fermi gas to date is T/Tp = 0.05 0j. This 
temperature may be obtained by adiabatic switching of a 
molecular BEC of temperature T = 0.5 Tbec with fi mo \ ~ 
ksT, which is routinely achieved for atomic BEC's. A 
temperature of T = 0.25 Tbec with /x mo i — feflT, as for 
example achieved for an atomic Bose condensed gas in 
Ref. O, would give T/T F = 5 x 10" 3 . 

This offers exciting possibilities for obtaining deeply 
degenerate Fermi regimes, such as observing the BCS 
transition in the weakly interacting regime kp\a\ <C 1. 
The corresponding critical condition is obtained by sub- 
stituting into E q. 11911 the BCS critical temperature cal- 
culated in Ref. |19j : 

T B cs/T F = aexp(-7r/2fc F |a|), (20) 

where a = e y ir~ 1 (2/e) 7 / 3 ~ 0.277, 7 being Euler's con- 
stant. Strictly speaking, in an isotropic harmonic poten- 
tial, Eq. (|20|l is only valid in the limit fcslBCS ^ ^ UJ - 



However, for kp\a\ = 1/2 or 1/3, it gives a result accu- 
rate to 30% for /cbTbcs — 27iw [23- In the right axis 
of Fig. n is shown the critical value (fci?|a|) c ensuring 
Tfinai = Tecs, as a function of T B ec/T: {k F \a\) c is on the 
order of one half, a regime already accessible by present 
experiments. It should be noted that our treatment is 
exactly valid only in the regime kp\a\ -C 1, as we have 
used BCS theory. The values of (kp \a\) c that we find 
are not extremely small compared to unity; we therefore 
expect corrections to both the equation of state of the 
normal phase and the critical temperature. Our result 
must therefore be considered as an estimate. 

Fermi gases have been predicted to be sensitive to heat- 
ing due to particle losses, much more so than Bose con- 
densates |23j . One may ask what the limit is on the acces- 
sible final temperature set by this heating during the adi- 
abatic tuning of a. Assume a loss rate 7i oss of atoms due 
mainly to collisions with the background gas, and let t be 
the duration of the adiabatic ramp. Using the derivation 
of Ref. , and assuming that 7loss t < 1 and (T/Tp) 2 < 
1, the temperature increase due to loss may be esti- 
mated as {T/Tp) i ncrcasc — a/tioss*/ (27r 2 ) . One may take 
t ~ l/[<y co \i(T/Tp) 2 ], where 7 co n = navp is the clas- 
sical collision rate involving the scattering cross-section 
a ~ Ana 2 . Here n is the total density, vf = Tikp/m is 
the Fermi velocity, and a factor of (T/Tp) 2 is included to 
account for the Pauli blocking for each of the two compo- 
nents ||. Then, supposing (T/Tp)n na _\ ~ 10~ 2 , one finds 
the condition 

[7 loss /(27r 2 7coll )] 1/4 < (T/T F ) final ~ 10~ 2 (21) 

in order to avoid a loss-induced increase in tempera- 
ture during the switching time. For the experiment of 
Ref. E3, 7i oss < 10~ 2 s"\ and Eq. CD is satisfied for 
a > 1000ao and 11 > 5 x 10 12 cm -3 , where ao is the Bohr 
radius. 

In order to satisfy the criterion of adiabaticity, ther- 
mal equilibrium must be maintained while tuning the 
scattering length from positive to negative values. This 
requires that the thermalization rate, i.e., the rate of 
elastic collisions, be much larger than the rate of change 
of a For a > 0, thermalization is due to colli- 

sions between non-condensed molecules. As the size of 
the non-condensed cloud is close to the size of the con- 
densate for ksT rsj jUmoii the effective collision rate is 
thus 7cff = 7SS 1 Nt/NjboI, where N T /N mo i is the non- 
condensed fraction and 7"^/ = n mo i<T mo \VT , with vt the 
thermal velocity and <7 mo i = 87ra 2 riol the scattering cross 
section between two molecules. For T/Trv.c, = 1/2, the 
non-condensed fraction is (T/Tbec) 3 = 1/8 [1J|. Assum- 
ing a ~ 1000ao and the mass of 6 Li, even for molecular 
densities n mo \ as small as 10 12 cm' 3 , this gives a mini- 
mum thermalization time of ~ 30 ms. For a < 0, the rate 
of thermalization is strongly reduced by Pauli blocking, 
and one has 7 cr f = j C oii(T/Tp) 2 0. At the critical value 
of kp\a\ ~ 1/2, and for a total density of 5 x 10 12 cm" 3 , 
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adiabaticity requires that the change of the scattering 
length occurs on a time t > 300 ms. This time is longer 
than the inverse oscillation frequency of lithium atoms in 
typical magnetic traps, so that the change in a will not 
induce macroscopic oscillations in the gas. 

It is thus far not possible to obtain a perfect balance 
of spin states in a Fermi degenerate gas 0. It may 
therefore be supposed that a small fraction of fermions 
will remain unpaired for positive scattering lengths, due 
to a lack of partners, and will coexist with the molecu- 
lar condensate. Fermions have an entropy proportional 
to T, as opposed to T 3 in the case of bosons. So even a 
small number of atoms could be expected to dominate the 
total entropy of the atom-molecule mixture at low tem- 
peratures, thereby interfering with the proposed cooling 
scheme. This problem can be avoided by simply remov- 
ing the excess fermions from the system for a > with a 
correctly tuned laser pulse, which requires that their ex- 
citation frequencies are sufficiently different from those 
of the molecules [22| . 

In conclusion, we have shown that adiabatic switch- 
ing of a molecular BEC allows one to obtain a deeply 
degenerate Fermi gas of temperatures on the order of 
T ~ 10~ 2 Tf- This suggests a way to achieve a BCS tran- 
sition in the weakly interacting regime. It is important to 
note that our thermodynamic approach does not require 
a knowledge of the equation of state of the system in the 
intermediate strong coupling regime where kp\a\ > 1. 
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Note added: Since submission of this paper, long- 
lived BEC's of weakly bound molecules have been ob- 
served [25| , and a strongly interacting fermionic conden- 
sate in the BEC-BCS crossover regime was created |26j . 
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